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In this work we investigate the dynamical properties of a mixture of mutually interacting spherical 
molecules of different masses and sizes. From an analysis of the microscopic laws governing the 
motion of the molecules we derive a set of non-local self-consistent equations for the singlet phase- 
space distribution functions. The theory is shown to reproduce the hydrodynamic equations for 
, the densities of each species, the total momentum and the local temperature. The non ideal gas 

^ ■ interaction term is separated into a contribution due to the repulsive part, which is treated by 

means of the revised Enskog theory for hard spheres, and an attractive contribution treated within 
] the random phase approximation. The present formulation accounts for the effects of the density 

D . and velocity inhomogeneities both on the thermodynamic and transport properties of the fluid. 

pL^ ' In a special limit, where one species is massive and diluted, the theory leads to a description which 

is formally identical to the dynamic density functional equation governing the time evolution of a 
colloidal system. The derivation also determines the dependence of the friction coefBcient, appearing 
in the dynamic density functional theory, on the microscopic parameters of the solvent. However, 
the predicted value takes into account only the coUisional contributions to the friction and not the 
f"|~| , Stokes friction of hydrodynamic origin, suggesting that velocity correlations should be incorporated 

Q ■ in a more complete treatment. 

?3 ■ I. INTRODUCTION 

' Spatially inhomogeneous systems on mesoscopic length scales can generate properties which do not appear in bulk 
(~| , materials, offering new perspectives for future applications. A variety of tools have been utilized to study the properties 
O ' of fluids near substrates and liquid interfaces in terms of molecular forces, ranging from new experimental techniques 
^ O ^ I to numerical algorithms and theoretical approaches [H-ll, S [l^l • 

It is an understatement to say that, since its appearance in 1979, the Manifesto 0] of density functional theory (DFT) 
has strongly influenced the studies in the field of classical inhomogeneous fluids. In this approach the equilibrium 
^ 1 density profile is determined by a functional derivative of a non-local Helmholtz free energy functional. An exact 
0^ ' theorem states that such an equilibrium profile minimizes the grand potential of the system and is unique. This fact 
, renders the method extremely appealing and provides a great help in finding good approximate Helmholtz functionals 
■ in a variety of cases. 

• I Regarding out of equilibrium systems, we do not have such useful theorems [8| and a time dependent extension 
of these methods, the dynamical density functional theory (DDFT), has been applied on phenomenological grounds 
Q : As shown by means of comparisons with numerical simulations, DDFT is able to capture the over-damped 

. Although some 
can be included 



^ dynamics of suspensions, but does not account for the richer dynamics of molecular fluids 11- 

• • ■ hydrodynamic aspects, such as the presence of a drift in the solvent or inertial corrections 

' ^^"^ DDFT, a full treatment of these effects requires a different approach. In order to capture the isothermal 

hydrodynamic behavior, describing the long length and long time scale behavior of fluids, one must consider the 
; I [ momentum density in addition to the mass density. Hydrodynamics, developed much earlier than Kinetic theory and 
Ci ' without any knowledge of the underlying microscopic structure of liquids, but using only phenomenological arguments 
and conservation laws flf\, represents a universal theory of fluid behavior. With this remark we want to stress the 
fact that, in constructing a dynamical description of a fluid valid both at molecular and macroscopic scales, we must 
be consistent with the assumptions which rendered classical hydrodynamics such a successful theory of liquids. In 
other words, symmetries and mass and momentum conservation laws must be exactly preserved even when necessary 
and unavoidable approximations are introduced. 

When a liquid is treated as a mechanical continuum, it is relatively simple to find a closed set of governing hydrody- 
namic equations for the relevant fields with the help of the so-called phenomenological constitutive relations. However, 
when the physical inhomogeneities become of the same order of magnitude as the molecular scales, one is forced to 
adopt a microscopic description [l^. The challenge is to extend macroscopic concepts such as pressure, viscosity, 
diffusivity and thermal conductivity to situations, where matter is confined to narrow spaces or is heterogeneous. All 
these properties assume a non-local dependence on the controlling fields and this fact renders the theory particularly 
challenging. 
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The present paper is organized as follows: in Sec. |TT] we introduce the model and the evolution equations for the 
distribution functions of the multi-component system. Numerical solutions of the transport equations are possible by 
means of the Lattice Boltzmann technique [l^ as discussed elsewhere (20l - [23l | . In Sec. IIIII we analyze the predictions 
of the theory concerning the equilibrium and non equilibrium behavior of the model and give the relevant formulae to 
compute the pressure and the transport coefficients. In Sec. IIVI we mimic a colloidal suspension by specializing the 
discussion to a hard-sphere binary mixture, composed of a low-concentration heavy species and a high concentration 
light species. From this non primitive model we obtain an equation for the evolution of the heavy species alone very 
similar to the DDFT equation. This heuristic derivation also shows the limits of approaches based on the neglect of 
velocity correlations. Finally in Sec. |V]we present some conclusions and perspectives. 



II. THEORY 



We consider an M-component fluid, whose species, denoted by the label a = 1, M , have masses and interact with 
pair additive, centrally symmetric potentials U°'^{r) and external forces Fq,. We describe the statistical evolution of 
the system by means of the so-called BBGKY hierarchy of equations [13], whose first M equations involve the singlet 
phase space distribution functions, /"(r,v,i) and the two particle distribution functions, f2^{T,T',w,w',t) : 

dtrir, V, + V V/"(r, V, t) + ■ (r, v, = V f^"^(r, v, t) (1) 

where the interaction term fl"^ is given by: 

f}"'^(r,v,t) = -Lv„-y' dV J dr'VrU"P{r-v')ff{r,r',^y,t). (2) 

The evolution equations for the singlet distributions involve the two particle distributions and these in turn involve 
the three-particle distributions. In order to make some progress one needs to approximate these higher order dis- 
tributions, using some physically motivated prescription. A popular closure ansatz allowing to reduce the BBGKY 
equations [23] to a closed set of equations for the singlet distributions is the following: 

/f (r, r', V, v', t) g^^ir, r', t)/" (r, v, t)f{r', v', t), (3) 

where gap(j^r\t) is the local equilibrium pair distribution function. The approximation (j3]) incorporates static 
correlations correctly, but neglects velocity correlations, through the lack of velocity dependence in (^^^(r, r', <). We 
shall comment later on the consequences of such an ansatz. We assume that the interaction potential between two 
particles can be separated into a short-range strongly repulsive and a longer range attractive contribution, so that the 
interaction process results in a combination of almost instantaneous, hard core collisions and small velocity changes 
induced by the weak attractive field: 

17"''(r,v,t) = f7,"4(r,v,t) +l]^f,(r,v,t). (4) 

The repulsive part of ^^ep is treated using the Revised Enskog theory of Ernst and van Beijeren [25l - l29| for hard-sphere 
mixtures of diameters CTq^, which neglects velocity correlations for two particles about to collide, as in Boltzmann 
theory, but includes the configurational correlations resulting from the finite size of the particles via the pair correlation 
function at contact: 

r!?fp(r,v",v'5,t) =a2^y" d^^^ J dse(s • v,^)(s • v„^) x 
[g^p{r, r - a^ps, t)r{r, v", t)/(r - a„^s, t) 

-g^p{v, r + a„0S, t)/" (r, v", ^/^(r + '^a^s, i)}, (5) 
where Q{x) is the Heaviside function, v^^ — (v" — v'^), while v" and are scattered velocities given by 



2m" 
to" -1~ mP 



v^=v^-t- , , , (s-v.^)s (6) 
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and s is the unit vector directed from particle a to particle /3. The quantities gapi^, r ± fc^s) are the inhomogeneous 
hard sphere pair correlation functions evaluated when the particles of species a and /3 are at contact distance CTq^ = 
(ctcq + cr^^)/2. The attractive term has the random phase approximation [23| expression: 

where G"'' are the molecular fields 

G"'^(r,t) = -J dr'n^(r',%„^(r,r',i)V.C/:,i(r-r') (8) 

and the number density of species f3. Even with approximations ([5]) and (O the full solution of eq. ([T]) is exceedingly 
difficult and can be found only in some special cases of limited interest. In order to encompass this problem many 
authors adopted simpler forms of the interaction term. Among these forms, a very popular recipe is represented by 
the BGK model [3d| consisting in replacing 51"^ by a relaxation term. For collisions between particles belonging to 
the same species one chooses: 

^%GK = -o^Urir, V, t) ~ ^"(r, V, t)) (9) 
and for collisions between unlike particles: 

^B%K = -u;oa(r(r,v,i) - ^"(r,v,t)) (10) 
where "0" is the local equilibrium distribution for species a: 

0"(r, V, i) = n"(r, ^)[^]^/^ exp( ..^^ ) • dD 

and the functions '0" and the frequencies have to be modeled to account for collisions between particles of different 
species fsTj . 

The BGK approximation is based on the idea that the system after few molecular collisions reaches a state of 
local thermodynamic equilibrium, where the distributions depend on space and time only through the hydrodynamic 
variables of the system, p"(r, t), u(r, t), T(r, t), partial mass densities, average velocity and temperature, respectively. 
These variables are defined in terms of the distributions as: 

p"(r,t) = m"n"(r,i) = m" / dv/"(r,v,i), (12) 



and 

= 5 ™" / '^^^^ " ^fr{r. V, t) (14) 

with total density given by 'n{r,t) = "-"(''i barycentric velocity by u(r, t) = p"(r, t)u"(r, t)/p(r, t) and 
global density p(r, t) — 0- With the help of the BGK ansatz one easily derives the hydrodynamic equations 

for the conserved variables /3"(r, t), pu(r, i) and the energy and via the Chapman- Enskog analysis [S^ the transport 
coefficients. Unfortunately, the method gives an oversimplified picture of the thermodynamic properties of the system, 
which turn out to be the same as those of an ideal gas, since f^^^^ does not contribute to the pressure or to the surface 
tension. A partial remedy to such a situation was introduced by Shan and Chen and other authors [33l - [35| . To take 
into account the contribution of the interactions to the equation of state they included a self-consistent term, named 
the pseudo-potential. This ad-hoc adjustment allowed to consider the hydrodynamic properties of non ideal gases 
by means of the so-called Lattice Boltzmann method (LBM) [l^. The pseudo-potential is a conservative force and 
does not determine the transport coefficients, which only depend on the adjustable parameters uJap- In the language 
of the present article the Shan-Chen method is a particular choice of G"^ with the peculiarity that it contains both 
attractive and repulsive contributions, an assumption which is somehow in conflict with the standard methods of 
liquid state theory, where one treats separately and on a different basis attractive and repulsive forces fs^ . 
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The full collision RET operator ([5]) being a non-linear functional of couples the different velocity moments of 
the distributions thus rendering analytic work very hard, unless one adopts a suitable truncation scheme. 

Any satisfactory approximation must retain the physical symmetries and conservation laws which are incorporated 
in the microscopic representation ([5|) and (O . To achieve that goal Dufty et al. [s^, HI] separated the contributions 
of il"^ to the hydrodynamic equations from those affecting the evolution of non-hydrodynamic modes, by projecting 
the collision term onto the hydrodynamic subspace spanned by the functions {l,v, u^} and onto the complementary 
kinetic subspace: 



f^"^ = Vhydron"^ + (/ - Vhydro)^''^ (15) 

with 



n,..o^^"^(r,v,0^ , \ ^"(r,v,0 ( (v-u).%(r,i) 

\ ^ 3kBT{r,t) ^'^ y^^'') 



(16) 



and 

/ \ , / 1 \ 

17"^(r,v,i). (17) 



C"'3(r,t) I - / dv f TO"(v-u) I 



The vanishing of the first element of the array in the l.h.s of eq. ([T7| expresses the conservation law of the number of 
particles of each species in a collision. 

Due to the splitting (jisp , the orthogonal part of the collision term will not appear explicitly in the balance equations 
for the hydrodynamic variables (see eqs. ([20| -([23 |) ). This fact suggests a simple approximation for (/ — Vhydro)^"'^ , 
where one replaces the exact expression by a BGK-like relaxation time term, having the property of vanishing under 
the application of the operator Vhydro- 

- Vhydro)n''^{r, V, t) ~ -co[r{r, V, t) - (r, v, t)], (18) 
where a; is a relaxation frequency and 

TO"(u"(r, t) ~ u(r, t)) ■ (v - u(r, t)) 



V'I(r,v,t) =V"(r,v,t){l + 



(19) 



The factor multiplying the Maxwellian in eq. ([T^ serves as to "orthogonalize" the term —Ld[f°' — 4>±] to the collisional 
terms proportional to C and B, so that the BGK contribution does not explicitly affect the balance equations. 

We construct, now, the hydrodynamic equations by projecting the equations ([T} onto the hydrodynamic subspace. 
Multiplying by {l,m"v,rn"(v — u)^/2}, integrating over velocity and summing over components, we obtain the 
following set of balance equations: 



5tp(r,i)+V - (^p(r,t)u(r,t)j=0, (20) 

which represents the continuity law of the mass density. In addition we have the continuity equations for each species: 

9tp"(r,t)+V - (p"(r,t)u"(r,i))=0. (21) 

We also find the momentum balance equation 

dt [p(r, t)uj (r, t)] + V, (p(r, t)u, (r, t)u, (r, <)]) 

= -V./^f ) (r, t) + Y^ (r, t) + Cf (r, t) (22) 

a a/3 
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and the balance equation for the focal temperature: 
^kgnir, t) (dt + u,{r, i)V,) T{r, t) 

= -/^f)(r,i)V,t.,(r,t)-V.gf\r,t)+i?(r,t)+^p"(r,i)^(<(r,t)-«.(r,t)), 

(23) 

with B = 'Ylap , where we have introduced the kinetic part of the pressure tensor: 

Plf\r,t)=Y.m'^ [ dv{vr-u,){vf-u,)r{r,v,t) (24) 
and the kinetic part of the heat flux vector 

qW^^m" /'dv(v-u)^^^l^^r(r,v,i). (25) 

a 

Eqs. and (|23l) assume their standard hydrodynamic form, if we insert the relation (see ref.js^l) between the 
coUisional moments and the divergence of the collisional transfer contribution to the pressure: 

V,PiPir,t)=-Y,Cfir,t) (26) 

a/3 

and the relation between the collisional transfer contribution to the heat flux and the pressure: 

V,gP(r,i) +4"" ^(r,t)V,u,(r,i) = -B{v,t). (27) 

Notice two facts: a) at equilibrium the term ()18p vanishes, b) out of equilibrium it determines a fast relaxation of 
the distributions towards the local values of u{r,t) an T{r,t). 

In order to obtain explicit expressions for C and B we make a further approximation. We perform the integrals 
appearing in eq. pT]) by replacing the true distribution functions f°'{r,v,t) by the Maxwellian distributions, corre- 
sponding to average density n"(r, t), local velocity u"(r, f) and temperature T(r, t). Many years ago the same kind of 
approximation was used by Longuet-Higgins and Pople |39j to predict the transport coefficients of hard core systems. 
As shown in ref. [12] a useful representation of the last term in eq. (l22t is obtained by the following decomposition: 

J2 C^'^ir, t) = n"(r, t) (F"'"-'^(r, t) + F"'''™f (r, t) + F"''''*'=°"^(r, t) + F"'^(r, t)) . (28) 

a 

The first term in the r.h.s. can be identified with the force acting on the a-particles at r due to the influence of all 
remaining particles in the system, the gradient of the potential of mean force: 

Y'^^^f{v,t) = -kBTY,^lf5 [ dssg^p{r,r + a^ps,t)npir + a^ps,t)+J2G'"^i^,t), (29) 

where the integral over s is over the surface of a unit sphere. The second term is the average drag force experienced 
by particles of species a when moving with velocity with respect to the remaining species having velocities u^: 



X 



(u"(r, t) - u'^(r, t)) ■ J dsssgapir, r + ctq/js, t)n'^{r + a^ps, t), 

(30) 

where ^Map is the reduced mass: ^lap — {m"m^)/{m" + m^). The third term is a viscous force due to the presence of 
gradients in the velocity field: 



(r,i) = ^2ai,^^^^^i^M>ll))x 

dssgai3{v,v + aai3S,t)n'^{v + aaps,t)s ■ (u'^(r + CTq^s) - u'^(r)). 



(31) 
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while the last term takes into account the contribution to the force due to temperature gradients: 



13 

n^{v + a^fiS,t)kB[T{v + (T^fiS,t)-T{v,t)]. (32) 



Finally, the energy integrals are given by the expression: 

B{Y,t) = kBT{Y,t)^a'ip j dsgafj{Y,Y + aa[is,t)n" {r,t)n'^ {y + aa[is,t) x 



s • [u^(r + a„^s, t) - u^{y, t)] - ^— • [u^'{y, t) - u(r, t)] 



to" 



.s.K(r,t)-u(r,t)] + 



2 2^i^pkBT{Y, t) [T(r + CT„0S, t) - T{y, t)] 



m°' + mP ^ ^ ' ' ^ ' " m^ + ml^M n T{Y,t) 

-n°'{Y,t)G°'^{Y,t) ■W°'{Y,t), 



(33) 



where w"(r,i) — u"(r, t) — u(r,t). Notice that eq. represents a decomposition of the total force in dissipative 
and non dissipative terms. The first ( eq. (1^^]) ) and the fourth term ( eq. (|3H)) are non dissipative and have 
opposite sign under time- reversal with respect to the partial momentum current n"u", that they induce via the 
balance equation eq. f (l^^ ). The dissipative terms eg. ([50)1 and ea. dSTI) . instead, have the same parity as the current 
(40| . Moreover, dissipative forces have an equilibrium counterpart, being related, as shown in ref.[23'|, to the intrinsic 
chemical potential of the individual components: 

F",-/(r,i)^„VM?„*(r,i), (34) 



while dissipative forces vanish at equilibrium. Eq. p4p represents a direct connection between the effective force 
featuring in the DDFT and the present approach and this fact opens the possibility of transferring to systems out 
of equilibrium the vast knowledge accumulated in the last twenty years concerning effective interactions in colloidal 
solutions and liquid mixtures at equilibrium (4l| - |43j . On the other hand, the dissipative forces cannot be derived 
from a functional derivative of the Helmholtz free energy functional, J-{{n"'}). In other words, the DDFT does not 
give information about the transport coefficients of the system, which are originated by non-equilibrium processes not 
accounted for by this approach. 

Once the self-consistent fields C and B have been specified, it is possible to solve eqs. ([IJ and study fluids under 
a variety of inhomogeneous conditions, using the LBM to achieve a numerical solution of eq. ([T|). We have shown 
in a recent series of papers that the main features of a structured fluid can be captured up to moderate packing 
fraction [2^ by a suitable extension of the LBM. In the present paper we shall not discuss the numerical aspects of 
the problem, but investigate further some theoretical issues. 



III. EQUILIBRIUM AND NON EQUILIBRIUM PROPERTIES OF BULK SYSTEMS 

In order to recover the thermodynamic properties we impose the conditions of global equilibrium. These require 
that all the velocities it" are equal and all hydrodynamic fields are time independent and, with the exception of the 
densities n", are spatially uniform. With T =constant u" =constant, eq. ([2^ . the momentum balance condition, 
reduces to the hydrostatic equilibrium condition. The r.h.s. of eq. (|26p relating the equilibrium part of the potential 
contributions to the pressure tensor can be separated into a repulsive, Pl^^Y,t) and an attractive term P°-p^{Y,t). 
The pressure can be obtained by applying the following formulae originally derived by Kirkwood [13, El] ■ 

Pl^^^{Y,t) = ^ M:,3^ I dsJ.^^. X 
1 

dXgapir - (1 - A)cra/3S, r + AcrQ/3s)7i"(r - (1 - X)aafss)n^{Y + Xdafss) 



(35) 
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and 



Q/3 



R OR 



f dXga^pir + (1 - A)R,r - AR)7i"(r + (1 - A)R)?i''(r - AR). 

"'0 



In the case of a uniform system the total bulk pressure, p''"'*^^ obtained by summing the kinetic and potential 
contributions and taking all the fields to be constant, is diagonal and isotropic: 

p,ui, ^ksTY^n'^{l + ^Y: -U^'9o.,M -^E^'J dRR^g^p{Rf-^^) . (36) 

In addition, the surface tension of a planar interface can calculated from Kirkwood and Buff formula 

/oo 
dz[PN{z) - Pt{z)], (37) 
-oo 

where the subscripts N and T indicate normal and tangential components of the pressure tensor, respectively. 

Non equilibrium properties. The kinetic coefficients can be considered to be the sum of two contributions: the first 
due to the instantaneous transmission of momentum and energy across the bodies of molecules upon collision and the 
second resulting from the distortion of the Maxwell-Boltzmann distribution induced by the presence of viscous and 
heat flows. Hereafter, we briefly report the calculation of the kinetic coefficients along lines similar to those discussed 
in ref. [2l|. 

Collisional contributions to transport coefficients. While kinetic transport, prevailing at low densities, is originated 
by the movement of the particles carrying a certain amount of momentum and energy, collisional transport, dominant 
at higher densities, is due to the transfer of momentum and energy from one particle to the other during collisions. 
Following steps similar to those of references [l^, [2l| , we derive the collisional shear viscosity and heat conductivity of 
the mixture. To simplify the derivation we assume that the densities and the temperature are uniform , ua^ ns — const 
and T(r, t) = Tq and u" = = u and that the velocity field is a shear slowly varying over distances of the order of 
the molecular diameter: 

u(r,t) = (0,wy(a;,0),0). (38) 
We compute the collisional contribution to the viscosity using the relation: 

a/3 

where we employed cq. (|26D for the first equahty and the constitutive relation 

p{C) ^ p(C) ^ (C) ^ ^1 UQ) 

for the second equality. With the help of eqs. ((28|) and (|3ip . after expanding to second order in (Tq,/?, we obtain: 

'l''^^ = 4 E V^^^i^^^^^^^Ugcf^n'^n^- (41) 

a/3 

Similarly, if we assume — constant , u" = u*^ and small temperature gradients, we can derive the thermal 
conductivity. We use the constitutive relation 

q(*=')(r,i) = -A('^)Vr(r,t), (42) 

compare with relation ([27| and expand B{r,t) ( formula ([33])) to second order in CTq^ with the result: 

A(^) = hsV^k^Y. y ABg.pn'^n^- (43) 
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The present results are consistent with the theory of transport coefficients put forward for mono-disperse hard-sphere 
systems by Longuet-Higgins and Pople more than half a century ago ^39| . It is based on the fact that in hard sphere 
systems even a local Maxwellian approximation to the distribution functions is able to account for the collisional 
transfer contribution to the transport coefficients. According to formulae (jH]) and it is clear that within the 
present approximation the mean field term, G"^, does not contribute to the transport coefficients. 

Kinetic contributions to transport coefficients. The self-consistent Longuet-Higgins and Pople gaussian approxima- 
tion employed above gives an expression for the transport coefficients at high densities, but does not give an expression 
for the kinetic contribution to these quantities. In order to obtain the full expression of the transport coefficients we 
apply the Chapman-Enskog analysis [32[ . The kinetic contribution to the transport coefficients can be derived upon 
neglecting the interaction terms C and B in the transport equation. The solution of the transport equation ^ can 
be written as: 

r(r, V, t) « V"(r, V, t) + Sr{r, v, t) (44) 
In this case eq.(IT]) to first order in Sf" reads: 

(9t + V . V + ■-^)r (r, V, t) = -uj6r{r, v, t) . (45) 

Substituting (|44p into (|45|) and taking the appropriate velocity moments we obtain the following Euler-like ( i.e. 
without dissipative effects) hydrodynamic equations for the mixture. These are: the continuity equation for each 
species 



dn" d{n°'uj) 



dt dxj 



= 0, (46) 



the global momentum conservation: 

^m"n"[^ . + ;^F.,(r)n"(r,i), (47) 

a ^ ^ a. 

and the temperature equation: 

^ (r, (f + ..|^) = - E t)T'^.^.,. (48) 



Using again the solution (|44)) in eq. (|45l) and eliminating the time derivatives of the hydrodynamic fields with the 
help of eqs. (P5 1) - P5)) we obtain the first order correction to the distribution function: 



;((.. - u.)iv, - .,) - -(V - + [m-'—^ - 2^BT)i-^-} 



A:rT' 



(49) 

where the first term contributes to the diffusive current, the second term to the viscous flow and the third to the heat 
flow. With the help of cq. (H^ . we compute the off diagonal part of the kinetic stress tensor (z ^ j): 



where the kinetic viscosity of the mixture is 



(50) 



Finally, we obtain the kinetic contribution to the heat flux 



q(^) ^ E / ^^(^ - u)i^— ^5/"(r, V, t) = -A(^) VT, (52) 
where the kinetic contribution to the heat conductivity is: 

2uj ^-^ m° 



aW-^V— . (53) 

9,,, m.a ^ ' 
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IV. HEURISTIC DERIVATION DDFT FROM HYDRODYNAMICS OF AN ASYMMETRIC BINARY 

MIXTURE 



Classical dynamic density functional theory (DDFT) is a widely used tool for studying the dynamics of suspensions 
of colloidal particles in a solvent [ZB - lioj . Within this approach only the density of the colloidal particles appears 
explicitly, whereas the solvent appears through the friction coefficient, or friction tensor. As a little remark, we wish 
to show that in the case of an isothermal hard sphere mixture in which the c (for colloidal) component is very diluted 
and the s (for solvent) particles are light im'^ /m'^ « 1) it is possible to derive an equation for the density of species 
c, which is a special instance of the DDFT equation. 

Being lighter, the s particles have higher thermal velocities and reach equilibrium faster. It is possible to integrate 
out their degrees of freedom and represent their influence by a viscous drag force, proportional to the velocity of the 
solute particles, plus a random force mimicking the effect of their random motions. 

From the continuity equation (|21l) for the density of c particles we have: 



ft/(r, t) + V • (u(r, t)p^{v, t))+V- ^!M^!M(u=(r, t) - u%r, t)) = 0, (54) 

where we separated the advection term from the diffusion term, and u'^ and are the velocity of the colloidal particles 
and of the solvent, respectively. We now assume that the inertial term in the momentum equation for the colloidal 
particles can be neglected with respect to the molecular forces and the motion is not accelerated. In other words, we 
consider a low Reynolds number regime and obtain: 

" -Vn'=(r,i) +r'=(r,t) +F^'"^(r,i) +F'='''™s(r,<) +F"'"''^(r,i) ~ (55) 



n=(r,t) 



where we used eq. (f28|) . If, in addition, we neglect the shear force, with the help of eqs. (|30|) and (|55|) we eliminate 
the drag force in favour of (u"^ — u''): 



u^(r,i)-u^(r,<)): 

3 TO* 



■ (/cbTV lnn'=(r, t) - F'='"^(r, t) - F"(r, i)) . (56) 



8p {2nfi,skBTy/^alg, 

We define the local chemical potential, /.t^(r,t), as the sum of the internal forces (see ref. [l!]): 

Vn^ir^t) = fcBTVlnn^(r,i) - F'''"'^r,t), 
substitute into the advection-diffusion equation ([5l)) . and take the diluted limit -nF- jrf << 1 and u ~ u'* 

dtTf(Y, t)+V ■ (u(r, t)n%r, t)) = • [n^(r, t) (Vfi^r, t) - F=(r)^ 

where the friction coefficient, 7, is given by the expression: 

1 3 



7 8n« ^TTm^kBTa^sdc 



(57) 



(58) 



We, now, observe that eg. dST]) is formally identical to a DDFT equation for the c species in a velocity field u(r,t). 
As — ?► the gradient of the chemical potential /i^ approaches the ideal gas value, kBTWn''{r,t), so that eq. ((57)) 
becomes a linear advection-diffusion equation for the field n^, with a diffusion coefficient given by: 

D='-^, (59) 

7 

to be interpreted as a fluctuation-dissipation relation between 7 and D . Apparently, such a result sounds correct, 
appealing and satisfactory, but contradicts Einstein's theory of Brownian motion stating that the diffusion coefficient 
for a large and massive solute molecule of diameter CTcc immersed in a solvent of much smaller and lighter molecules 
is related to the solvent viscosity by the Stokes-Einstein relation: 

DsE = (60) 

hydro 
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Einstein's theory relies on the hypothesis that a diffusing spherical body moves among solute particles as a macroscopic 
sphere does in a viscous incompressible continuum fluid. One first needs to solve the Stokes equation for the flow 
of the fluid around the sphere, with the condition of no slip at its surface. From the flow one computes the total 
stress acting on the sphere surface and finally the drag force, which results proportional to the (hydrodynamic) radius, 
crcc/2, of the object and to the viscosity, r/s, of the solvent according to the Stokes formula: 

pcMrag ^ _^^^^^^(u- _ = _3^7?.aec(u^ - u). (61) 

The Enskog relation (|58|) and the Stokes-Einstein relation (|6T|) display different scalings with respect to the physical 
properties of the c particles. The relation 

J hydro — 3TTT]sCrcc 

for the friction coefficient shows a linear dependence on the solute diameter CTcc, but no dependence on its mass, while 
the Enskog approximation considerably underestimates the friction coefficient, which depends on the reduced mass of 
the solute-solvent pair and has a quadratic dependence on the solute diameter. Only for very heavy solutes the mass 
dependence of Enskog friction coefficient vanishes (see eq. ([5S)) ) . The numerically computed friction coefficient shows 
a crossover to the Stokes-Einstein result only for large values of the mass ratio at fixed size ratio. A full merging 
between the two approaches, in spite of repeated attempts, has not been achieved so far ^-^]. 

Clearly, the Enskog microscopic picture, where the interactions occur via collisions with the solvent particles, and 
the Einstein picture, which treats the solvent molecules as a continuum and the effect of the solute-solvent interactions 
by means of the no slip boundary condition, account for different mechanisms. Einstein's theory is correct for large 
solute particles, where the typical ratios between the radii of suspended colloidal particles and solvent molecules 
range in the interval 30 — 3000 and the mass ratio is in the interval 10'' — 10^". These large differences justify the 
the idea of complete separation of time scale between the motion of the solute, which occurs in a time of th e order o f 
Tc = m'^/(37r77*(Tcc), and the motion of the solvent characterized by the Enskog collision time ~ n'^ al^yj UbT / , 
the first being 5 orders of magnitude larger [s^. However, for small colloids the coUisional and hydrodynamic regimes 
arc not well separated. A crossover should occur when the size and the mass of the solute become comparable to 
those of the solvent molecules. In this case eq. (pT|) is not expected to be valid. The short time direct collisions 
between the tracer particle and the solvent particles are indeed taken into account by Enskog's theory, which, however, 
misses the long time behavior determined by correlated re-collisions of the tracer with the same solvent molecule (H^j ■ 
The velocity auto-correlation function (VACF) of the c particles in Enskog's theory displays an exponential decay 
which is valid in the short time region, but does not show the algebraic inverse power law tails characterizing the 
long time behavior of the VACF. The pure exponential behavior is the result of the molecular-chaos assumption, 
i.e. of a Markovian representation of the dynamics, well describing the large mean free path regime. The RET 
is able to capture the formation of positional correlations, but not the build-up of velocity correlations among the 
particles, at the origin of the Stokes-Einstein behavior of the diffusion coefficient. It is our opinion that methods, 
which truncate the BBGKY hierarchy at the level of one-particle distribution functions, are not apt to account for 
the renormalization of the transport coefficients due to the presence of velocity correlations, the ultimate cause of 
hydrodynamic interactions. With respect to this problem mode coupling theory (MCT) is able to predict the correct 
behavior by including at the same time binary collision processes, coupling with transverse current of the solvent and 
with density fluctuations jss} . 

Finally, it is worth to mention that a simpler way to bridge Einstein with Enskog theory as the size or the colloids is 
reduced to that of the solvent particles has been described in ref. [24;, 5(| , where the motion of the colloid is described 
by a Langevin equation with a friction coefficient which is non-local in time in order to reproduce a non Markovian 
behavior. Although, such an approach is frankly phenomenological it is able to predict the backscattering effect. 

V. CONCLUSIONS 

In this paper we have considered how microscopic methods, utilized to study the non-equilibrium evolution of 
colloidal solutions and well accounting for their configurational properties, can be extended to investigate the dynamics 
of liquid mixtures. We have adopted a strategy based upon kinetic methods and shown that the resulting equations 
correctly describe the hydrodynamic behavior of the density, momentum and energy transport. We have found 
similarities with the DDFT equations, and this not surprising since the common basic ingredient of both methods is 
the free energy which in one way or another features in the theory. In DDFT we directly use the free energy functional, 
and in the kinetic approach we employ its functional derivative, the inhomogeneous positional pair correlation function, 
g{r, r'). Finally, to link further the two descriptions we have introduced an heuristic derivation of the DDFT equation 
from the kinetic equations for a hard-sphere mixture. Our result seems to indicate the possibility of extending this 
procedure to mixture of higher physical interest. 



11 



Acknowledgments 

It is a pleasure to dedicate this paper to Bob Evans, whose work has deeply influenced the development of condensed 
matter in the last forty years. I had the privilege to work with him in the H.H. Wills Physics Laboratory at the 
University of Bristol. 



T.M. Squires and S.R. Quake, Rev.Mod.Phys. 77, 977 (2005). 
R.B. Schoch, J.Han, P. Renaud, Rev.Mod.Phys. 80, 839 (2008). 

W. Sparreboom, A. van den Berg and J. C. T. Eijkel, Nature Nanotechnology 4, 713, (2009). 

R. Evans, in Fundamentals of Inhomogeneous Fluids, ed. D. Henderson, Dekker, New York, (1992), ch. 3. 

M. Rauscher and S. Dietrich, Annual Review of Materials Research 38, 143 (2008). 

H. Bruus, Theoretical Microfluidics, Oxford University Pres., New York, 2008. 

R. Evans, Adv. Phys. 28, 143 (1979). 

G. K. L. Chan and R. Finken, Phys. Rev. Lett. 94, 183001 (2005). 
A.J. Archer and R. Evans, J. Chem. Phys. 121, 4246 (2004). 

U. Marini Bettolo Marconi and P.Tarazona, J. Chem. Phys. 110, 8032 (1999) and J.Phys.: Condens. Matter 12, 413 

(2000) . 

A. J. Archer, J. Phys.: Condens. Matter 17, 1405 (2005). 

A. J. Archer, J. Chem. Phys. 130, 014509 (2009). 

M. Rauscher, J. Phys.: Condens. Matter 22, 364109 (2010). 

F. Cecconi F, F.Diotallevi, U. Marini Bettolo Marconi and A. Puglisi, J. Chem. Phys. 120 35 (2004). 

U. Marini Bettolo Marconi and P. Tarazona, J. Chem. Phys. 124, 164901 (2006); U. Marini Bettolo Marconi, P.Tarazona 
and F. Cecconi, J. Chem. Phys. 126, 164904 (2006). 

U. Marini Bettolo Marconi and S.Melchionna, J. Chem. Phys. 126, 184109 (2007). 

G. K. Batchelor An Introduction to Fluid Dynamics, Cambridge University Press, Cambridge (1967). 
K. Huang, Statistical Mechanics John Wiley, New York, (1963). 

S. Succi, The lattice Boltzmann equation for fluid dynamics and beyond, 1th edition , Oxford University Press, Oxford, 

(2001) . 

S. Melchionna and U. Marini Bettolo Marconi, Europhys.Lett 81, 34001 (2008). 

U. Marini Bettolo Marconi and S.Melchionna, J. Chem. Phys. 131, 014105 (2009). 

U. Marini Bettolo Marconi and S.Melchionna, J. Phys.: Condens. Matter 22, 364110 (2010). 

U. Marini Bettolo Marconi and S.Melchionna. larXiv : 1009.5600 (submitted to J.Chem.Phys) (2010). 

J. P. Hansen and I.R. McDonald, Theory of Simple Liquids Academic Press, Oxford, (1990). 

H. van Beijeren and M.H. Ernst, Physica A 68, 437 (1973) and 70, 225 (1973). 

M. Lopez de Haro, E.G.D. Cohen and J.M.Kincaid, J.Chem.Phys. 78, 2746 (1983). 

J. Karkheck and G. Stell Phys. Rev. A 25, 3302 (1982). 

J. M. Kincaid, S. Prez, and E. G. Cohen Phys. Rev. A 38, 3628 (1988). 

J. Karkheck, E.Martina and G. Stell Phys. Rev. A 25, 3328 (1982). 

P. L. Bhatnagar, E. P. Gross, and M. Krook, Phys. Rev. 94, 511 (1954). 

V. Sofonea and R.F. Sekerka, Physica A 299, 494 (2001). 

S. Chapman and T.G. Cowling The mathematical Theory of Non-uniform Gases 3rd ed. Cambridge University Press, 
Cambridge 1970). 

X. Shan and H. Chen, Phys. Rev. E 47, 1815 (1993); X. Shan and H. Chen, Phys. Rev. E 49, 2941 (1994). 
X. Y. He and L. S. Luo, Phys Rev E 56, 6811 (1997). 

S. Chen and G. D. Doolen, Annual Review of Fluid Mechanics 30, 329 (1998). 
J. D. Weeks, D. Chandler and H. C. Andersen, J. Chem. Phys. 54, 5237 (1971). 

A. Santos, J.M. Montanero, J.W. Dufty and J.J. Brey Phys.Rev. E , 57, 1644 (1998). 
J. W. Dufty A. Santos, and J. Brey Phys. Rev. Lett. 77, 1270 (1996). 

H.C. Longuet-Higgins and J. A. Pople, J.Chem.Phys. 25, 884 (1956). 

P. M. Chaikin and T. C. Lubensky, Principles of Condensed Matter Physics, Cambridge University Press, Cambridge 
(1995). 

T. Biben and J. P. Hansen, Phys. Rev. Lett. 66, 2215 (1991). 

R. Roth and R. Evans and A. A. Louis, Phys. Rev. E 64, 051202 (2001). 

B. Gtzelmann, R. Evans, and S. Dietrich, Phys. Rev. E 57, 6785 (1998) 

J. A. McLennan Introduction To Nonequilibrium Statistical Mechanics, Prentice Hall, New York (1988). 

J. Z. Wu and Z. D. Li, Annual Review of Physical Chemistry 58, 85 (2007). 

J.G. Anero and P.Espanol, Europhys.Lett 78, 50005 (2007). 

M. Rauscher, J.Phys.: Condens. Matter 22, 364109 (2010). 

M. Rex, H. Lowen, and C. N. Likos, Phys. Rev. E 72, 021404 (2005). 

P. Tarazona and U. Marini Bettolo Marconi, J. Chem. Phys. 128, 164704 (2008). 



12 



[50] J.T. Hyncs, R. Kapral and M. Weinberg, J. Chem. Phys. 70, 1456 (1979). 

[51] J. R. Schmidt and J.L. Skinner, J. Chem. Phys. 119, 8062 (2003). 

[52] Wokyung Sung and G.Stell, J. Chem. Phys. 80, 3350 (1984) and J. Chem. Phys. 77, 4636 (1982). 

[53] L. Bocquet, J. Piasecki and J.P Hansen, J. Stat. Phys. 505, 76, (1994). 

[54] S. Musharaf Ali, A. Samanta and S. K. Ghosh, Chem. Phys.Lett . 217, 357 (2002). 

[55] S. Bhattacharyya and B. Bagchi, J. Chem. Phys. 106, 1757 (1997). 

[56] B.J. Berne, J.P. Boon and S.A. Rice, J. Chem. Phys. 45, 1086 (1966). 



